Entropy exchange and entanglement in the Jaynes-Cummings model 
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The Jaynes-Cummings model is the simplest fully quantum model that describes the interaction 
between light and matter. We extend a previous analysis by Phoenix and Knight (S. J. D. Phoenix, 
P. L. Knight, Annals of Physics 186, 381). of the JCM by considering mixed states of both the 
light and matter. We present examples of qualitatively different entropic correlations. In particular, 
we explore the regime of entropy exchange between light and matter, i.e. where the rate of change 
of the two are anti-correlated. This behavior contrasts with the case of pure light-matter states 
in which the rate of change of the two entropies are positively correlated and in fact identical. 
We give an analytical derivation of the anti-correlation phenomenon and discuss the regime of its 
validity. Finally, we show a strong correlation between the region of the Bloch sphere characterized 
by entropy exchange and that characterized by minimal entanglement as measured by the negative 
eigenvalues of the partially transposed density matrix. 



Quantum entropy was first formulated by von Neu- 
mann [lj as an extension of the Gibbs entropy in classical 
statistical mechanics. The foundations of modern infor- 
mation theory were established by Shannon |2| using a 
definition of entropy which is essentially identical to the 
Gibbs entropy. Similarly, quantum entropy plays a cen- 
tral role in the theory of quantum information Q Q 
Modern definitions of quantum entropy include those of 
Renyi [| and Tsallis 0j. 

In recent years, entropy of quantum systems has 
been discussed in the context of entanglement. The 
Horodecki's derived Renyi entropy inequalities and 
discussed their implications for non-locality. They 
showed that the entropic inequalities are violated by any 
pure entangled state. They also argued that by consider- 
ing Renyi entropy inequalities one can limit teleportation 
U°f several states. Cerf et at gave an interest- 
ing interpretation to independent, correlated, and entan- 
gled two qubit systems by considering conditional and 
mutual entropies. For example, a maximally entangled 
(non-separable) two qubit system possesses negative con- 
ditional entropies and excessive mutual entropy (which is 
related to unaccessible information), while a maximally 
correlated (separable) two qubit system possesses zero 
conditional entropies and maximal classical mutual en- 
tropy. Phoenix and Knight [TTj | explored the dynamics 
of a single cavity mode in resonance with a single atom 
by comparing von Neumann entropy, Shannon entropy, 
and photon number variances. 

In this paper we focus not on the entropy of the indi- 
vidual subsystems, but on the entropy correlations be- 
tween subsystems in a bipartite system. This question 
has been almost completely overlooked in the literature, 
and for a good reason: for a bipartite system that starts 
in a pure state, the partial entropies of the two subsys- 
tems are equal at all times. This result is a consequence 
of the equality between the eigenvalues of the partial den- 
sity matrices in a pure bipartite system, which in turn is 
a consequence of the Schmidt decomposition. In this case 
the correlations are trivial. However, there are two qual- 
ifications in obtaining these trivial correlations. First, 



that the system is indeed bipartite. Second, that the joint 
state of the two subsystems is pure. If either of these con- 
ditions is not met the correlation between the entropies 
of the subsystems is nontrivial. In particular, in this pa- 
per we focus on the question of whether and under what 
conditions it is possible to get anti-correlated behavior 
of the entropy of the subsystems. This is interesting be- 
cause it implies that there is an entropy transfer process 
going on, consistent with classical thermodynamic con- 
cepts but opposed to the default result obtained for en- 
tangled bipartite systems. In this paper, for dcfinitencss, 
we focus on entropy correlations between an atom and 
a single quantized cavity mode in the framework of the 
JCM. We present qualitatively different entropy corre- 
lations between the atom and the field, and we demon- 
strate a regime of entropy exchange between them, both 
numerically and analytically. 

This paper is arranged in the following manner. Sec- 
tion I is a brief introduction. It is devoted to a descrip- 
tion of the density matrix of bipartite systems, the defi- 
nition of quantum entropy, and the issue of separability 
and entanglement of bipartite systems. In this section 
we also review the JCM Hamiltonian in the context of 
quantum entropy. Section II is devoted to calculating 
entropy correlations for mixed states of light and matter 
with qualitatively different behaviors. We give examples 
of entropy exchange between light and matter, and ex- 
plore the parameter range of this regime. In Section III 
we give analytical proof of the entropy exchange effect 
and an analysis of the regime of its validity. In Section 
IV we apply two different entanglement tests to deter- 
mine whether the atomic-field system is entangled for all 
types of entropy correlations. Section V concludes. 



I. INTRODUCTION 

A bipartite system is described by a density matrix of 
a C m ® C n Hilbert space. The partial density matrix of 
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one part is obtained by tracing over the other: 

PA(B) =Tr B ( A ){p AB ). (1) 

The entropy of a quantum system is given by the von 
Neumann entropy 0: 

S = -k B Tr{p\np). (2) 

The purity of a quantum system is given by Tr{p 2 ), and 
it is bounded: < Tr{p 2 ) < 1. Purity is related to the 
q = 2 Tsallis entropy 0: 

S 2 = l-Tr(p 2 ). (3) 

Qualitatively, purity oppositely tracks the von Neumann 
entropy: an increase in the von Neumann entropy is par- 
allel to a decrease in purity (or an increase in the Tsallis 
entropy) . 

A bipartite quantum system is considered separable if 
it can be written as [T^ ]: 

PAB = ^P i p i A®P B , (4) 

where P % > 0, P l — 1; an d P A m) are individual par- 
tial density matrices. A system that can not be factored 
into the form above is said to be entangled. In practice, 
determining whether a general C m ® C™ bipartite system 
can be factored into the form above is very hard. Thus, 
several tests have been introduced in order to determine 
whether a system is entangled. One test originates from 
quantum information theory, and it relies on calculating 
conditional entropies. Conditional entropy indicates the 
entropy of one subsystem after measuring the other, and 
it is given by: 

S(A\B) = S AB - S B . (5) 

A bipartite system with at least one degree of freedom 
having negative conditional entropy is entangled. There- 
fore a necessary condition for separability is that the con- 
ditional entropies are positive. Alternatively, one can 
calculate the mutual entropy. Mutual entropy indicates 
the entropy shared between the two subsystems, and it 
is given by: 

S(A:B)=S A + S B -S AB . (6) 

Mutual entropy is bounded: 

0<S(A:B)<2min[S A ,S B ]. (7) 

A bipartite system whose mutual entropy is excessive: 
min[S A , Sb] < S(A : B) < 2min[SA, S B ] is entangled. 
Therefore a necessary condition for separability is < 
S(A : B) < min[S A , Sb}- 

A second powerful test introduced, originally by Peres, 
[l3| relies on partial transposition. Partial transposition 
is a blockwise transposition of a matrix and it is given 
by: 

Pia,j0 = Pip,]*- ( 8 ) 



A system whose partially transposed density matrix is 
negative is entangled. Therefore a necessary condition for 
separability is the positivity of the partially transposed 
density matrix (the Horodecki's showed that it is also 
a sufficient condition for separability of C m £g> C" systems 
with ran < 6). Kraus and coworkers [l5j analyzed sys- 
tems supported on C 2 <£> C N . They gave necessary and 
sufficient conditions for separability of systems with pos- 
itive partially transposed (PPT) density matrices. They 
showed that if the rank of density matrix is equal to N, 
the density matrix is separable. However, other crite- 
ria they devised are very complicated to implement for a 
generic density matrix. 

The JCM is the simplest fully quantum model that 
describes the interaction between light and matter. The 
model consists of a single quantized two level atom in- 
teracting with a single quantized electromagnetic cavity 
mode under the rotating wave approximation (RWA) and 
the dipole approximation. The resonant JCM Hamilto- 
nian |lq | in the interaction representation is given by: 

H = h\((j + a + (T-a ] ). (9) 

For a system that starts in a pure state partial entropies 
of the field and the atom are equal at all times for a 
system that starts in a pure state As already indi- 
cated above, this result is a consequence of the equality 
between the eigenvalues of the partial density matrices 
in a bipartite system, which in turn is a consequence of 
the Schmidt decomposition. Since the partial entropies 
fluctuate together in time, their sum does not conserve 
the total entropy. 

S a +Sf>S af . (10) 

Thus, in this case partial entropies are not an additive 
(extensive) quantity. Nevertheless, knowledge of partial 
entropies gives information about the dynamics of the to- 
tal system. The question that we address here is whether 
we can find states of light and matter for which the par- 
tial entropies are quasi-additive. Mixed states of light 
and matter can be natural candidates in this respect. 



II. RESULTS 

We consider various combinations of initial pure and 
mixed states of light and matter. The mixed state for the 
electromagnetic field is chosen to be a Planck distribution 
for a single cavity mode: p f (Q) = Y^=o p n\n) {n\, P n = 
(71+")"+ 1 ' wnere ^ is the average number of photons in the 
cavity. The mixed state for the atom has the following 
general form: p a (0) = P e |e)(e| + (1 - P e )\g) (g\, < P E < 
1. When the cavity mode is in exact resonance with the 
atomic transition, one can derive an analytic expression 
for the propagator and obtain analytic expressions for the 
full and partial density matrices at any time. The general 
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solution for the full density matrix has the following form: 



Paf(t) = e-^ Paf (0)e 



AHt 



(11) 



where H is the interaction Hamiltonian mentioned above. 



By expanding U(t) 



-iHt 



in a Taylor series one obtains 



the analytic form of the propagator: 



U(t) = 



■ sinf XtV aat ) 

-I — v , — '-a 

V aaT 



- sin( Xt V af a) f / \ . / + \ 

-i — v ; , — L a' cos(XtVa'a) 

V a< a 



To simulate the evolution of the entropic quantities, we 
truncated the set of Fock states that compose the Planck 
distribution at some nt, where Y^nLo Pn — 1- The ac- 
curacy of the results we obtained was tested by adding 
more Fock states to the Planck distribution to see if the 
simulated values changed. 

Entropy change 
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FIG. 1: Partial entropy change plots AS = S - 5(0) (field 
dotted line, atomic solid line). Top: atom initially in the 
excited state and a weakly excited thermal field, p(0) = 
(|e)(e|) a ® (^2n=o Pri\n){n\\ n — 0.1)/. Bottom: atom ini- 
tially in the ground state and a weakly excited thermal field , 
P(0) = (lff)<sl)a®(£~ -P„|ri>(n|) / ;n = 0.1) / . The dash-dot 
curve represents the sum of atomic and field partial entropy 
changes, and it is quasi-conserved. 

In Fig. ^the changes in atomic and field partial en- 
tropies (S(t) — 5(0)) vs. the dimensionless quantity Xt 
are plotted. In Fig. la the atom is initially excited and 
the field is in a weakly excited thermal state (n = 1). 
This case is similar to a situation where both the field 
and atom are initially in a pure state, because both par- 
tial entropies rise and lower together (although they are 
not equal). In Fig. lb the atom is initially in the ground 
state and the field is in a weakly excited thermal state. 
It is clear that there is entropy exchange between the 
atom and field (although the exchange is not complete). 
The sum of the field and atomic entropy changes (green 



curve) is quasi-conserved as seen by its amplitude of fluc- 
tuations, which is significantly smaller than each of the 
partial entropy changes. 

When the field is excited more, i.e. n — 1, or n — 
10, no substantial entropy exchange between the atom 
and field occurs. In particular, when the field is highly 
excited (n > 10) and the atom is close to a pure state, 
there is a sharp and rapid collapse of the atomic purity 
(sharp and rapid rise in the von Neumann entropy) with 
no substantial revival. 

In order to determine which initial atomic states can 
exchange entropy with a weakly excited thermal field we 
considered different initial atomic states by discretizing 
the Bloch sphere with a longitudinal angle (— J < 9 < ■§ , 
for positive 6 the atom is more excited) and a Bloch vec- 
tor length (0 < r < 1), disregarding the azimuthal angle 
((f)). We introduce a time averaged entropy exchange pa- 
rameter defined by: 



.AS, 



«(/)• 



•A5 



/(«) 



(12) 



where A5 a (/) = S a tf)(tj)— 5 a (/)(tj-_i), and in the numer- 
ator the smaller of the two instantaneous partial entropy 
changes was substituted. The time averaged entropy ex- 
change parameter is suitable for quantifying the extent 
of entropy flow between a weakly excited thermal field 
and an atom due to the oscillatory nature of the partial 
entropies. Furthermore, the entropy exchange parameter 
is bounded: — 1 < P < 1; when it tends to —1 there is 
high degree of entropy exchange, and when it tends to 
+1 the two partial entropies rise and lower together. 

In Fig. [5] the entropy exchange parameter is plot- 
ted for a weakly excited thermal field and various ini- 
tial atomic states. The dark blue region in parameter 
space is characterized by significant entropy exchange 
(P < —0.8); and in this region the atom is mostly in 
the ground state (ground state probability is more than 
80%). Thus entropy exchange occurs most effectively 
when the atom is initially close to the ground state. 
An example for an almost complete entropy exchange 
is given in Fig. |3 where the changes in atomic and field 
partial entropies are shown for an atom initially in the 
state p a (r = 0.7,9 = -f) (P « 0.99). In this case 
the sum of the two partial entropies is almost completely 
conserved as seen by the amplitude of fluctuation of the 
sum of the partial entropy changes, which is two orders 
of magnitude smaller than each of the partial entropy 
changes. 



III. ANALYTICAL DERIVATION OF THE 
ENTROPY EXCHANGE EFFECT 



We now provide an analytical derivation for the en- 
tropy exchange effect. We first note that contours of 
significant entropy exchange (P < —0.8, Fig. |2J) center 
around a point on the central vertical axis of the Bloch 
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Entropy exchange parameter 
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FIG. 2: Entropy exchange parameter for a weakly excited 
thermal field and various initial atomic density matrices, char- 
acterized by the longitudinal angle 9 and the Bloch vector 
length r. The dark region in parameter space is character- 
ized by significant entropy exchange (P < —0.8). Note that 
this region is centered around a 'fixed' point where the partial 
density matrices are stationary (black arrow; see text). 
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FIG. 3: Partial entropy change (field dotted line, atomic 
solid line, atomic+field dash dot curve) for a weakly ex- 
cited thermal field (n = 0.1) and atom initially in the state 

p a (r = \J~^\ & = "§)• Entropy exchange between the atom 
and field is almost complete (P ~ 0.99). 

By substituting eauation ll^l into eauation ll4l one obtains: 



sphere. This is a 'fixed' point where the partial density 
matrices are stationary, and it corresponds to a situation 
where the field temperature is identical to the atomic 
Boltzman temperature as we shall now show. Consider 
a situation where: 



Pe _ Pn+1 
~P ~ ~~P 



— P k R T 



n + 1 



(13) 



where to is the transition frequency. The full density 
matrix where both the atom and field are in mixed states 
(with no initial coherence) is given by: 



Paf(t) = U I Q 





Pa ®Pn 



where (P n ) mn = S m nPn- The partial density matrices 
are in diagonal form: 



Pa{t) 



(A{t) 
^ B(t) 

PT® = C(t) + D(t), 



(14) 



where 



A(t) = P e p n cos 2 Ki) + Pg Yl P «+i sin 2 (/3„+it), 

B(t) = J P e ^P n _ 1 sin 2 K_ 1 i)+P g ^P„cos 2 ( / 9 n i), 

C(t) = F e (P„cos 2 (a I ^)+F n _isin 2 K_ 1 i)), 

D(t) - P g (P n+1 sin 2 (/3 n+1 t) + P n cos 2 ((3 n t)), 



a n -i = Pn = Xy/ri. 



(15) 



Pa{t) 



Pa(0) - 



Pe 
P g 

P n - 



(16) 



For example for a weakly excited thermal field with (n) = 
0.1, this point is located at r = |; 8 = — -| (as indicated 
by the black arrow in Fig. 

We now proceed to analyze the vicinity around the 
fixed point, where substantial entropy exchange occurs. 
For a weakly excited thermal field the atomic fixed point 
is close to the ground state (see Fig. |2J|. Therefore, we 
consider an initial state where an atom in the ground 
state interacts with a thermal cavity mode: /Q Q /(0) = 

(\9)(9\)a® (J2n=o p n\ n )( n \)f- In tnis case the full density 
matrix is given by: 



n ,m - f*W0)sl -iS lPf (o)cl 



where C\ = cos ( Xt VaJa), and Si 
partial density matrices are in diagonal form: 



siniXtvaaA mi 

— ^= — a. the 



Pa(t) = 



(Y. P n+i sin 2 '(Pn+lt) 



J^Pn COS 2 {(3 n t) 



p] n (t) = P n cos 2 (p n t) + P n+1 sm 2 ((3 n+1 t). 



(17) 



Since the purity is qualitatively similar to the von Neu- 
mann entropy we can analyze the time dependence of 
the partial purities. The time derivatives of the partial 
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purities are given by: 
d(Tr(pl)) 



n =2J2 Pn + 1 sin 2 (/?n+it) p n+iPn+i sin(2/3 n+ iO 

- 2 Pn COS 2 {/3 n t) PnPn M^M 

d(Tr(p 2 f )) „ 

^ = 2^(P n+1 /? n+1 sin(2/?„ +1 t) - P n /? n sin(2/? n i)) 

• (F„cos 2 (/3„i) + -P«+isin 2 (/3 n+1 i)). (18) 

One can rearrange the time derivatives in equation 1181 
in terms of PiPj multiplying sums of oscillatory func- 
tions. For a weakly excited thermal cavity (n = 0.1) the 
PiPj decrease by approximately an order of magnitude 
for every increment in cither i or j. Approximating the 
derivatives of the partial purities by the first non vanish- 
ing term (PqPi = 0.0751) one obtains: 



vanishing term (Pq 
dt 



d{Tr{p})) 



dt 



0.8264) one obtains: 
— P 2 oi$ sin(4o!o£) 
— Pq ccq sin(4aoO- 



(22) 



The above approximation shows that in this case the par- 
tial purities oscillate together in time. The leading terms 
in the time derivatives of the partial purities are identical. 

When the cavity mode is excited thermally with more 
photons, i.e. n > 1, the terms in the time derivatives 
of the partial purities rearranged according to PiPj no 
longer decrease by an order of magnitude as i, j increase. 
Now many terms contribute substantially to the partial 
entropy behavior; since these terms are not identical in 
both degrees of freedom, observing entropy exchange be- 
comes increasingly less probable as n increases. 



d(Tr(pl)) 
dt 

d(Tr(p})) 
dt 



-2P Pi/3! sin(2/3it) 
2P P 1 /3 1 ewppit). 



(19) 



The above approximation shows explicitly that there is 
entropy exchange between a weakly excited thermal cav- 
ity mode and an atom initially in the ground state. The 
leading terms in the time derivatives of the partial pu- 
rities have the same functional dependence on time but 
with opposite signs. 

Consider now an initial state of an excited atom inter- 
acting with a thermal cavity mode: p a f(0) = (|e)(e|)„ (8 
(S^Lo Pn\ n )( n \) f- I n this case the full density matrix is 
given by: 



IV. ENTROPY CORRELATIONS AND 
ENTANGLEMENT 

We consider now the relationship between entropy 
correlations developed between the atom and field dis- 
cussed in previous sections, and entanglement of the 
various atomic-field states. One would expect that as 
Sa + Sb —> Sab the system becomes separable. This is 
motivated by the following formula from statistical me- 
chanics Il7l: 



Sab — Sa + Sb + k) 



r AB 



(23) 



Paf(t) 



C 2Pf (0)Cl iC 2Pf (0)Sl 
-iS^pfCl S 2 p f (0)Sl 



where C2 = cos( XtVaaJ), and So = 2Hii^*2Ji) a t_ Again 

V a< a 



the partial density matrices are in diagonal form: 

Pa(t) 
P f{t) 

The time derivatives of the partial purities are given by 



£P„cos 2 Ki) 

Y, p n-\ sin 2 (a n -it) / 

P n cos 2 (a n t) + P„_ x sin 2 (a„_ 1 t). (20) 



If Sab = Sa + Sb, then P^ B — P\P B (since the sum 
over the natural logarithm terms has to be zero). This 
means that the two subsystems are independent (if Sa + 
Sb — Sab then the two subsystems are weakly classically 
correlated). In regions of the Bloch sphere where the 
entropy exchange is almost complete (AS a —— —AS/): 



S a + S f ~ 5„(0) + 5/(0) = S o/ (0) = S af . 



(24) 



d(Tr{pl)) 
dt 



d{Tr{p))) 
dt 



-2 y P n a n sin(2a„i) } P n cos 2 (a n t) 



The last two equalities hold since the combined atomic- 
field system begins to evolve from a separable state, and 
since the evolution is unitary, respectively. Therefore, we 
would expect that in regions of the Bloch sphere where 
the entropy exchange is almost complete, there should be 
minimal entanglement. We now proceed with comparing 
2^^ Pn-i sin 2 (ct!„_i£) P n -xa n -x sin(2a n _ii)our measure for entropy exchange with various entangle- 
ment measures. 

As discussed in the Introduction, bipartite systems 
with excessive mutual entropy are entangled. We intro- 
duce an entanglement parameter R based on the ratio 
between mutual entropy and partial entropies: 



= 2 2_^(P„_ia„_i sin(2a„_i<) - P n a n sin(2a„i)) 



(P„ cos 2 (a„t) + P„_x sin 2 (a„_it)). 



(21), 



Rearranging the time derivatives in eauation l21l in terms 
of PiPj-, as in the case when the atom was initially in the 
ground state, and approximating them by the first non 



R 



S(a : /) 
min[S a ,S f ] 



(25) 
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Entanglement parameter, R=S(a:f)/min[S ,S ] 

1 a f 
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FIG. 4: Time averaged entanglement parameter based on the 
ratio between mutual entropy and partial entropies for dif- 
ferent initial atomic states coupled to a weakly excited ther- 
mal field (n = 0.1). The region in parameter space where 
R < 1 (dotted line) encompasses the entire region of entropy 
exchange (P < —0.8), but in addition includes significant area 
where there is no entropy exchange (even regions with P > 0; 
compare with fig. [3J 



FIG. 5: Entanglement parameter based on the negativity of 
the eigenvalues of the partially transposed density matrix for 
a weakly excited thermal field (n w 0.1) and various atomic 
states. The region of atomic parameter space with substan- 
tial entropy exchange (P < —0.8, bounded by a solid white 
line) falls within the region with minimal entanglement pa- 
rameter (dark) , which is dominated by anti-correlated partial 
entropies. 



It follows from equation [7| that < R < 2, with 
1 < R < 2 being a sufficient condition for entanglement 
and conversely < R < 1 being a necessary condition for 
separability. In fig. 0]the time averaged entanglement pa- 
rameter, R, is plotted for different initial atomic states 
coupled to a weakly excited thermal field. The region in 
parameter space with R < 1 (bounded by a dashed line) 
indicates that the necessary condition for separability is 
fulfilled for most of the evolution time. This region is 
large, occupying about 75% of the interior of the Bloch 
sphere. It encompasses the entire region of entropy ex- 
change (P < —0.8), but in addition includes significant 
area where there is no entropy exchange (even regions 
with P > 0). This is consistent with the observation 
that < R < 1 is a necessary but not a sufficient condi- 
tion for separability. Although the region defined by this 
necessary condition is too broad to be informative, note 
that in cases where there is substantial entropy exchange 
as seen in fig. [3 R approaches zero in fig. 0] 

We consider now the PPT test, where the existence of 
even a single negative eigenvalue of the partially trans- 
posed density matrix is a sufficient criterion for entangle- 
ment. Since the PPT test can only be applied for finite 
dimensional density matrices, in order to apply it to the 
Jaynes-Cummings model we need to truncate the infinite 
set of Fock states. We do this according to the following 



procedure: 

oo 

Pf (0) = £P„|n)<n| (26) 

n=0 

n t n f 

^J2 p n\n)(n\ + (1 -£p n )|n /+ i)<n /+ i|;(27) 

n=0 n=0 

in words, all residual probability from the truncation is 
placed in the nj + 1st state. The criterion for truncat- 
ing the infinite set of Fock states is that the partial en- 
tropies do not change by more than 10~ 14 . We find that 
no matter what the initial state of the atom is, there is 
always at least one negative eigenvalue of the partially 
transposed density matrix. However, further investiga- 
tion shows that this result is in essence due to the trunca- 
tion of the infinite Fock basis to a finite size. Specifically, 
we find that in all regions of the Bloch sphere there is 
always one negative eigenvalue of extremely small mag- 
nitude whose size is on the order of the probabilities tail 
of the infinite Fock basis that gets truncated. Since our 
primary interest here is in the original and not in the 
truncated Jaynes-Cummings model, we are inclined to 
ignore this negative eigenvalue of extremely small mag- 
nitude and to interpret the test of partial transposition 
in terms of the remaining negative eigenvalues. The re- 
maining negative eigenvalues are all of sizable magnitude 
and robust with respect to the size of the Fock basis. 

With that introduction, we define as the largest 
(in absolute value) negative eigenvalue of the partially 
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transposed density matrix, p 2 . We then introduce as an 
entanglement measure 

E = log(|A£|), (28) 

where the overbar signifies the time average. In figure 
|5]we plot the entanglement parameter, E, for different 
initial states for nf — 13. Clearly, there are two well- 
defined regions in atomic parameter space. The dark 
region corresponds to a single negative eigenvalue of p T2 
on the order of 10 _ 16 — 10~ 18 . This region is charac- 
terized by marginal entanglement, corresponding in all 
likelihood to completely separable evolution in the un- 
truncated Jaynes-Cummings model. Notice that the re- 
gion of atomic parameter space with substantial entropy 
exchange (P < —0.8 in fig. shown as a solid white line 
in fig. falls precisely within the region with minimal 
entanglement parameter, E. The lighter shade region in 
fig. [5] is characterized by two or more negative eigen- 
values of p T2 . We believe that this region corresponds 
to significant entanglement in the original, untruncated 
Jaynes-Cummings model. This region is clearly orthogo- 
nal to that of the region of entropy exchange. 

V. CONCLUSION 

We have explored entropy correlations between a quan- 
tized cavity mode and a single atom in the framework of 
the JCM by considering both pure and mixed atomic and 
field states. In particular, we explored the regime of en- 
tropy exchange between light and matter. We presented 
two qualitatively different entropy correlations. The first 
type of correlation is a case where both the atomic and 
field partial entropies fluctuate together in time. This is 
reminiscent of the case where both the atom and field 
start out in a pure state, and consequently their partial 
entropies are identical at all times. The second type of 
correlation is a case where the atomic and field partial 
entropies are anti-correlated. This implies that there is 
entropy exchange between the atom and the field. Since 
substantial entropy exchange occurs when the field is in 
a weakly excited thermal state (n = 0.1) we introduced 
an entropy exchange parameter in order to determine 
which initial atomic states can efficiently exchange en- 
tropy with the field. Substantial entropy exchange oc- 
curs when the atom is initially close to the ground state. 
We showed that contours of substantial entropy exchange 
center around a stationary point in the Bloch sphere. 
This point corresponds to a situation where the initial 



field and atomic (Boltzman) temperature exactly match. 
By analyzing the partial purities we derived an analytic 
approximation for the entropy exchange phenomenon. 
The change in the atomic and field partial purities have 
the same functional dependence on time but with oppo- 
site signs. 

It is natural to ask if there is any connection between 
the entropy exchange that we observe in this paper and 
a thermodynamic heat exchange process. Although we 
cannot rule out this possibility, we believe that the con- 
nection is unlikely. In classical thermodynamics, heat 
exchange is related to the transfer of energy that is in 
some sense degraded. In our case, because the entropy 
exchange is periodic and the system is small, there is 
no reason to believe a priori that the energy that is ex- 
changed is unrecoverable for pure work. Nevertheless, 
if one examines the ratio ggr^ , where dE/dt and dS/dt 
are the instantaneous rates of energy and entropy change 
of each subsystem respectively, there is some overall cor- 
relation with the initial temperature T of the field, in 
agreement with the classical result ^ = T. 

A well-established measure of entanglement is the test 
of negative eigenvalues of the partial transpose of the den- 
sity matrix (PPT test). Unfortunately, the PPT test can- 
not be applied directly to the Jaynes-Cummings model 
with its formally infinite Hilbert space for the field. We 
therefore applied the PPT test to a model in which the 
infinite Fock basis was truncated. We found that there 
is always at least one negative eigenvalue for the par- 
tially transposed density matrix. However, further inves- 
tigation showed that this result is due to the truncation 
of the Fock basis. Since our primary interest here is in 
the original and not in the truncated Jaynes-Cummings 
model, we essentially discard this negative eigenvalue of 
extremely small magnitude and interpret the test of par- 
tial transposition in terms of only the remaining negative 
eigenvalues. We find that the region of the Bloch sphere 
in which there are are no additional negative eigenval- 
ues maps very closely onto the region where where there 
is substantial entropy exchange between atom and field. 
This result is intuitively appealing, showing a strong neg- 
ative correlation between entropy exchange and entangle- 
ment. 
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